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^ ■ Abstract 



Vector-boson fusion processes are an important tool for the study of electroweak 
symmetry breaking at hadron colliders, since they allow to distinguish a light Higgs 
boson scenario from strong weak boson scattering. We here consider the channels 
WW — > ZZ and ZZ -^ ZZ as part of electroweak Z boson pair production in associ- 
ation with two tagging jets. We present the calculation of the NLO QCD corrections 
to the cross sections for pp^e"^e~ /x+/i~ -|- 2 jets and pp—>e'^e~u^i>^ + 2 jets via 
vector-boson fusion at order a^ a®, which is performed in the form a NLO parton-level 
Monte Carlo program. The corrections to the integrated cross sections are found to 
be modest, while the shapes of some kinematical distributions change appreciably at 
NLO. Residual scale uncertainties typically are at the few percent level. 



1 Introduction 

One of the primary goals of the CERN Large Hadron CoUider (LHC) is the discovery 
of the Higgs boson and a thorough investigation of the mechanism of electroweak (EW) 
symmetry breaking 012]. In this context, vector-boson fusion (VBF) processes have emerged 
as a particularly interesting class of processes. Higgs boson production in VBF, i.e. the 
reaction qq -^ qqH, where the Higgs decay products are detected in association with two 
tagging jets, offers a promising discovery channel [3] and, once its existence has been verified, 
will help to constrain the couplings of the Higgs boson to gauge bosons and fermions j3] . 

In order to distinguish possible signatures of strong weak-boson scattering from those of a 
light Higgs boson, a good understanding of WW -^ ZZ and ZZ — > ZZ scattering processes, 
which are part of the VBF reaction qq -^ qqZZ, is needed. This requires the computation 
of next-to-leading order (NLO) QCD corrections to the qq -^ qq ZZ cross section, including 
the leptonic decays of the Z bosons. Experimentally, very clean signatures are expected 
from the ZZ —>■ £~^i~i'~^i'~ decays in VBF with four charged leptons in the final state, the 
disadvantage of this channel being a rather small Z —>■ e+e" oi Z ^ /^^/^^ branching ratio of 
about 3%. The Z Z ^>- d,'^ d,^ vv channel, with two undetected neutrinos, on the other hand, 
results in a larger number of events due to the larger Z ^*vv branching ratio j3] . 

LO results for EW ZZ jj production in VBF have been available for more than two 
decades. The first calculations [H] were performed employing the effective W approxima- 
tion [7], where the vector bosons radiated off the scattering quarks are treated as on-shell 
particles and, therefore, kinematical distributions characterizing the tagging jets cannot be 
predicted reliably. In the following years, exact calculations for qq — > qq ZZ have been com- 
pleted, first without Z boson decay jS] , and then including leptonic decays of the Z bosons 
within the narrow width approximation [H]. 

We go beyond these approximations and develop a fully-flexible parton level Monte Carlo 
program, which allows for the calculation of cross sections and kinematical distributions for 
EW ZZ 23 production via VBF at NLO QCD accuracy. The program is structured in com- 
plete analogy to the respective code for EW W^W~ jj production presented in Ref. [TU] . 
Here, we calculate the t-channel weak-boson exchange contributions to the full matrix ele- 
ments for processes like qq ^ qq e^ e~ ^~^ fi~ and qq^qqe~^e~ p^p^ at 0{a^as). We consider 
all resonant and non-resonant contributions giving rise to a four charged-lepton and a two 
charged-lepton plus two neutrino final state, respectively. Contributions from weak-boson 
exchange in the s-channel are strongly suppressed in the phase-space regions where VBF 
can be observed experimentally and therefore disregarded throughout. We do not specifi- 
cally require the leptons and neutrinos to stem from a genuine VBF-like production process, 
but also include diagrams where one or two of the Z bosons are emitted from either quark 
line. Diagrams, where the final state leptons stem from a 7 -^ £"^£~ decay or non-resonant 



production modes, are also taken into account. Finite- width effects are fully considered. 
For simplicity, we nonetheless refer to the qq^qqi^i^i'^i'~ and qq -^ qq l^ (!.~ vv processes 
computed this way generically as "EW ZZ jp production. 

The outline of the paper is as follows. In Sec. |21 we briefly summarize the calculation 
of the LO and NLO matrix elements for EW ZZ jj production making use of the helicity 
techniques of Ref. JI]. Section El deals with phenomenological applications of the parton- 
level Monte Carlo program which we have developed. Conclusions are given in Sec. |3] 



2 Elements of the calculation 

The calculation of NLO QCD corrections to EW ZZ jj production closely resembles 
our earlier work for EW W^W~ production in association with two jets JU]- The main 
differences lie in the electroweak aspects of the processes, while the QCD structure of the 
NLO corrections is very similar. The techniques developed in Ref. jT^ can therefore be 
adapted readily and only need a brief recollection here. For simplicity, we focus on the 
e^e^ yU^yU^ decay channel in the following. The application of the basic features discussed 
for this case to the e^e~ v^v^ leptonic final state is then straightforward. 

The Feynman graphs contributing to pp -^ e^e^ fi^fi^ jj can be grouped in six topologies, 
respectively, for the 579 t-channel neutral-current (NC) and the 241 charged-current (CC) 
exchange diagrams which appear at tree level. These groups are sketched in Fig. Q for 
the specific NC subprocess uc—*uce^e~fi^fi^. The first two of these correspond to the 
emission of two external vector bosons V from the same (a) or different (b) quark lines. 
The remaining topologies are characterized by the vector-boson sub-amplitudes Lyy, Ty, 
Tyy and Tyy^, whlch describe the tree-level amplitudes for the processes VV -^ e^e~ /i"*"/^", 
V"— >e"'"e~ yU"'"/i~, VV ^* ^^ ii~ and VV ^>-e^e~. In each case, V stands for a virtual 7 or Z 
boson, and a and [3 are the tensor indices carried by these vector bosons. The propagator 
factors l/(g^ — 'fn'v +i'mvTv) are included in the definitions of the sub-amplitudes, which we 
call "leptonic tensors" in the following. Graphs for CC processes such as us ^dce~^e~ fi'^fj,~ 
are obtained by replacing the t-channel 7 or Z bosons in Fig. Q with W bosons. They 
give rise to the new lepton tensors L^+^_, T^+^_ ^ and T^+^_ for the sub-amplitudes 
W-^W- -^ e+e- fi-^fi", W+W- -^ e+e" and W+W ^ ^i+ ^i- . 

Contributions from anti-quark initiated t-channel processes such as uc^uce^e^ ^^^^, 
which emerge from crossing the above processes, are fully taken into account. On the other 
hand, s-channel exchange diagrams, where all vector bosons are time-like, contain vector- 
boson production with subsequent decay of one of the bosons into a pair of jets. These 
contributions can be safely neglected in the phase-space region where VBF can be ob- 
served experimentally, with widely-separated quark jets of large invariant mass. In the 
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Figure 1: The six Feynman- graph topologies contributing to the Born process 
uc ^ uc e'^ e~ n^ H^ . Diagrams analogous to (a), (d), (e) and (f), with vector-boson 
emission off the lower quark line, are not shown. 



same way, u-channel exchange diagrams are obtained by the interchange of identical final- 
state (anti) quarks. Their interference with the t-channel diagrams is strongly suppressed for 
tj^ical VBF cuts and therefore completely neglected in our calculation. 

For the treatment of finite-width effects in massive vector-boson propagators we resort 
to a modified version of the complex-mass scheme ^2); which has already been employed 
in Refs. P^ llOj . We globally replace vector-boson masses my with my — imyTv, without 
changing the real value of sin^ 9w This procedure respects electromagnetic gauge invariance. 
The amplitudes for all NC and CC subprocesses are calculated and squared separately for 
each combination of external quark and lepton helicities. To save computer time, only the 
summation over the various quark helicities is done explicitely, while the four distinct lepton 
helicity states are considered by means of a random summation procedure. 

The computation of NLO corrections is performed in complete analogy to Ref . PU] ■ For 
the real-emission contributions we encounter 2892 diagrams for the NC and 1236 for the CC 
processes, which are evaluated using the amplitude techniques of Ref. [TT] and the leptonic 
tensors introduced above. Singularities in the soft and coUinear regions of phase space are 
regularized in the dimensional-reduction scheme ^1] with space-time dimension d = 4 — 2e. 
The cancellation of these divergences with the respective poles from the virtual contributions 
is performed by introducing the counter terms of the dipole subtraction method [^. Since 
the color and fiavor structure of our processes are the same as for Higgs boson production 
in VBF, the analytical form of subtraction terms and finite collinear pieces is identical 
to the ones given in Ref. ^Hl- The finite parts of the virtual contributions are evaluated 
by Passarino-Veltman tensor reduction J7j, which is implemented numerically. Here, the 
fast and stable computation of pentagon tensor integrals is a major issue, which is tackled 
by making use of Ward identities and mapping a large fraction of the pentagon diagrams 
onto box- type contributions with the methods developed in JU]. The residual pentagon 
contributions amount only to about 1 %o of the cross sections presented below. 

The results obtained for the Born amplitude, the real emission and the virtual corrections 
have been tested extensively. For the tree-level amplitude, we have performed a compari- 
son to the fully automatically generated results provided by MadGraph y^, and we found 
agreement at the 10~^^ level. In the same way, the real emission contributions have been 
checked. For the latter, also QCD gauge invariance has been tested, which turned out to be 
fulfilled within the numerical accuracy of the program. The numerical stability of the finite 
parts of the pentagon contributions is monitored by checking numerically that they satisfy 
electroweak Ward identities with a relative error less than 6 = 1.0. This criterion is violated 
by about 3% of the generated events. The contributions from these phase-space points to 
the finite parts of the pentagon diagrams are disregarded and the remaining pentagon parts 
are corrected by a global factor for this loss. In Ref. [lOj we found that this procedure gives 
a stable result for the pentagon contributions when varying the accuracy parameter 6. 



3 Numerical results 

The cross-section contributions discussed in the previous section are implemented in a 
fully-flexible parton-level Monte Carlo program for EW ZZji production at NLO QCD 
accuracy, very similar to the programs for -ffjj, V 33 and W~^W~jj production in VBF 
described in Refs. jTH], jUj and jTU]. 

Throughout our calculation, fermion masses are set to zero and external b- and t-quark 
contributions are neglected. However, the code does allow the inclusion of 6-quark initiated 
sub-processes for NC exchange where internal top-quark propagators do not occur. For the 
Cabibbo-Kobayashi-Maskawa matrix, Vckm, a diagonal form equal to the identity matrix 
has been used, which yields the same results as a calculation employing the exact Vckm 
when the summation over all quark flavors is performed. 

We use the CTEQ6M parton distributions with as{mz) = 0.118 at NLO, and the 
CTEQ6L1 set at LO [19j. We chose mz = 91.188 GeV, mw = 80.419 GeV and Gp = 1.166 x 
10^^/GeV^ as electroweak input parameters. Thereof, uqed = 1/132.54 and sin^^iy = 
0.22217 are computed via LO electroweak relations. To reconstruct jets from flnal-state 
partons, the /ct algorithm [2011^ is used with resolution parameter D = 0.8. 



Partonic cross sections are calculated for events with at least two hard jets, which are 
required to have 

PTi>20GeV, \yj\<4:.5. (3.1) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four-vector sum of massless partons of pseudorapidity \ri\ < 5. The two reconstructed jets 
of highest transverse momentum are called "tagging jets". At LO, they are the flnal-state 
quarks which are characteristic of vector-boson fusion processes. Backgrounds to VBF are 
signiflcantly reduced by requiring a large rapidity separation of the two tagging jets. We 
therefore impose the cut 

^Vn = IVn - yj2\ > "^ ■ (3.2) 

Furthermore, the tagging jets are imposed to reside in opposite detector hemispheres, 

Vji ■ yj2 < , (3.3) 

with an invariant mass 

Mjj > 600 GeV . (3.4) 

These cuts render the LO differential cross section for ZZjj flnite, since they enforce 
finite scattering angles for the two quark jets. For the NLO contributions, initial-state 
singularities, due to coUinear q—^gq and g -^ qq splitting, are factorized into the respective 
quark and gluon distribution functions of the proton. Additional divergences, stemming 
from the t-channel exchange of low-virtuality photons in real-emission diagrams, are avoided 



by imposing a cut on the virtuality of the photon, Q^min ~ ^ GeV^. Events that do not 
pass this cut give rise to a colhnear q—'-qj singularity, which is considered to be part of the 
QCD corrections to p'j -^ ZZjj and not calculated here. 

In the discussion of the pp -^ i^i^i'^i'^ jj channel, we focus on the leptonic final state 
e^e~ fJ^^fJ^~ throughout. Results for the four lepton final state with any combination of elec- 
trons and/or muons (i.e. e~^ e~ fi~^ fi~ , e~^e~e~^e~ and fi'^ fi~ fi~^ fi") are obtained by multiplying 
the respective numbers for e~^e~fi~^fi~ by a factor of two. This procedure neglects very small 
Pauli-interference effects for identical charged leptons, however. Similarly, the e~^e~ v^v^ 
combination on which we concentrate, is related to an arbitrary two lepton plus two neu- 
trino final state in the case oi pp^i~^i~h'u jj production: Here, a factor of six is needed to 
take into account all neutrino species and two families of charged leptons. To ensure that 
the charged leptons are well observable, we impose the lepton cuts 

Pre > 20 GeV , \rie\ < 2.5 , ARje > 0.4 , 

mei > 15 GeV , ARu > 0.2 , (3.5) 

where ARjc and ARu denote the jet-lepton and lepton-lepton separation in the rapidity- 
azimuthal angle plane. In addition, the charged leptons are required to fall between the 
rapidities of the two tagging jets 

yj,min ^ "Ul ^ yj,max ■ V'^-'-'j 

In order to compute the full cross section for EW ZZjj production, contributions from 
the Higgs boson resonance, qq — »• Hqq —>■ ZZqq, as well as from the ZZ continuum have to be 
considered. A representative for the latter, which effectively starts at the Z-pair threshold, is 
obtained (for a Higgs mass below the Z-pair threshold) by imposing a cut on the four-lepton 
invariant mass of 

Mzz = V{Pi++Pi-+Pi'++Pi"y >mH + 10 GeV , (3.7) 

and correspondingly for the £"*"£" i/z/ final state. Since the contribution from the Higgs boson 
resonance has already been computed in Ref. ^^, we focus on the ZZ continuum in the 
following, if not stated otherwise, and assume mn = 120 GeV. 

The continuum cross section cTcuts for EW e^e~ u^j^u^jj production, within the cuts of 
Eqs. (I3.1|l - ()3.7j) . is shown in Fig. El The figure illustrates the dependence of the NLO 
cross section on the renormalization and factorization scales, /i/j and fip, which are taken as 
multiples of the Z mass, 

fJ'R = ^Rmz , fJ^F = ^F mz . (3.8) 



The LO cross section only depends on fip = ^p mz- By varying the scale factor ^i;' = ^ in the 

,LO 

cuts 



range 0.1 -r- 10, the value of a^^g changes by almost a factor of two, indicating the theoretical 
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Figure 2: Scale dependence of the total EW e^e v^yfiij cross section at LO and 
NLO within the cuts of Eqs. if^) -f3r7p for pp collisions at the LHC. The NLO 
curves show a^^^ as functions of the scale parameter ^ for three different 



cases: 



tJ-F = fJ-R = i'niz (solid red), fip = S,mz and fin = mz (dot-dashed blue), ^p = "mz 
and jjLR = ^mz (dashed green). The LO cross section depends only on fip (dotted 
black line). 



uncertainty of the LO calculation. The strong scale dependence is reduced substantially after 
the inclusion of NLO corrections. For 0"^^? "^^ study three different cases: ^p = ^b. = C, (solid 
red line), ^f = ^, ^r = ^ (dot-dashed blue line), and ^f = 1, ^_r = ^ (dashed green line). The 
latter curve illustrates clearly the weak dependence of cr^'jjf on the renormalization scale, 
which can be understood from the fact that as{^iR) enters only at NLO, while the LO cross 
section is completely independent of iir. Also the factorization-scale dependence of the full 
cross section is rather low, such that the variation of cr^^° with the scale parameter amounts 
to less than 2% for all cases in the interesting range 0.5 < ^ < 2. In the following, we fix the 
scales to /^f = /^ij = "^z- 

As a representative for the observables characterizing the tagging jets, we show the invari- 
ant mass distribution of the tagging-jet pair, da/dM,^^, for EW e~^e~ n^fi~ jj production in 
Fig. 121 The shape of the distribution at NLO differs from the respective LO result. This is 
emphasized in panel (b) of the figure, where we show the dynamical K factor, defined as 
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Figure 3: Invariant-mass distribution of the tagging jets in EW e'^ e^ fj,^ fi^ j j 
production at the LHC. Panel (a) shows the NLO (solid red) and the LO results 
(dashed black). Panel (b) displays the K factor as defined in Eq. i9..9j) . 



Due to the extra parton emerging in the real-emission contributions to the NLO cross section, 
the quarks which constitute the tagging jets tend to have smaller transverse momenta than 
at LO, thereby giving rise to lower values of their invariant mass. The transverse- momentum 
distributions of the tagging jets, per se, exhibit K factors in the range 0.8-i-1.4. On the other 
hand, angular distributions, such as the azimuthal angle between the tagging jets, display 
rather uniform K factors for the scale choice /i^ = /^_r = ^^z • For this particular scale choice, 
the total cross section is barely affected by the inclusion of NLO corrections, leading to a 
K factor of 0.99. 

Finally, we show the distribution of the invariant mass of the e~^e^ ^^fi^ system, which 
is given by Eq. (j3.7|) and can be fully reconstructed experimentally. It is very sensitive to a 
light Higgs boson, showing a pronounced resonance behavior for rriH ^ 800 GeV. For values 
of niH of the order of 1 TeV, the peak is diluted due to the large corresponding width of 
the Higgs boson {Th ~ 500 GeV) and the signal is distributed over a wide range in Mzz- 
Figure E] illustrates the resonance behavior of Mzz'- panel (a) shows the Mzz distribution for 
the continuum production of four charged leptons. Panel (b) displays the same observable, 
but now including the resonance contribution from a Higgs boson with rnn = 500 GeV. 

One remarkable feature of Fig. |3]is that LO and NLO results are virtually indistinguish- 
able, for the scale choice fip = fJ'R = ""^z- For the Mzz distribution, the invariant mass of the 
Z pair is another possible scale choice, which, however, would lead to substantially reduced 
LO differential cross section predictions at high values of Mzz- One finds, for example, a 
decrease by a factor of ~ 1.8 in da^'~'/dMzz at Mzz = 1-5 TeV when changing ftp = mz to 
Up = Mzz, which largely is due to an underestimate of the LO parton distributions at large 
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Figure 4: Panel (a): distribution of the four lepton invariant mass in EW 
e^e~ fJ-^f^^ jj continuum production at the LHC, within the cuts of Eqs. icl.l\} - ^[^71\ ) 
with niH = 120 GeV. Panel (b) shows the same observable when the contribution 
from a Higgs boson of mass mn = 500 GeV is included. In each case, NLO (red 
solid) and LO (black dashed) results are shown. 



Feynman-x. The NLO prediction, on the other hand, decreases by about 13% compared to 
our default choice fJ,F = f^R = ""^z, demonstrating the precision gained by including the NLO 
corrections. 



4 Conclusions 

In this paper we have presented first results for EW ZZ jj production at NLO QCD 
accuracy, obtained with a fully-flexible parton-level Monte Carlo program that allows for 
a straightforward implementation of realistic experimental cuts. The integrated cross sec- 
tions for the two processes pp -^ e+e~ /x"'"/i~ jj and pp -^ e~^e~ v^p^ jj were found to exhibit 
K factors around 0.99, which shows that higher-order corrections are under excellent control. 
Larger NLO contributions are obtained for some kinematical distributions with dynamical 
K factors in the range 0.8 ^ 1.4. These results hold for a default scale choice of /x = niz- 
Leading order results can change substantially, by up to a factor of 2, for other scale choices, 
while NLO results are very stable, demonstrating the value of the NLO corrections. An esti- 
mate of the theoretical uncertainty of the NLO calculation is provided by the scale variation 
of cross sections, within VBF cuts. It amounts to about 2% for integrated cross sections and, 
in extreme cases, up to 10% for distributions, when changing scales by a factor of two. Sim- 
ilar uncertainties are induced by the present errors on parton distribution functions, which, 
in the analogous case of Higgs boson production in VBF, were found to be about ±3.5% jT^ . 
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